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Abstract 

This report concerns the energy of a zero-temperature many-body system of spin i fermions interacting 
via a two-body potential with a free space infinite scattering length and zero effective range; the Unitary 
limit. Given the corresponding phase-shift S(k) = tv/2 a one-term separable potential is obtained by inverse 
scattering assuming a momentum cut-off A such that 5(k) = for k > A. The effective interaction in the 
many-body system is calculated in a pp-ladder approximation with Pauli-blocking but neglecting mean-field 
(dispersion) corrections; effective mass m* = 1. Using only the zero relative momentum component of this 
interaction the total energy is £ = 4/9 (in units of the fermigas), a result reported by several previous authors. 
Integrating the momentum dependent interaction over the Fermi sea this energy is revised to £ = 0.24. This 
result is independent of density and of the cut-off A if A >~ 3fc/. 

With m* 7^ 1 there is however a strong dependence on this cut-off. 

Including hh-ladders estimates give £ = 0.4 «-» 0.6, but a reliable result would in this case require a 
Green's function calculation. 

1 Introduction 

The properties of a dilute fermigas with large scattering length is of considerable theoretical as well as ex- 
perimental interest. Taking advantage of Feshbach resonances it is possible to magnetically tune the atomic 
scattering lengths, e.g.p] [2 [3]. Increasing the scattering length of fermions from — to +oo resulting in bound 
boson systems to explore the crossover from BCS to BEC has been reported by several groups. 

A theoretically related problem proposed by George Bertsch [3] is that of the energy of a dilute system of 
spin 1/2 fermions interacting via a zero-range, infinite scattering length interaction referred to as the Unitary 
limit. For such a system one would expect the existence of a constant £ being a function only of fundamental 
constants such that the total energy E = ^Efg where Efg is the uncorrelated Fermi-gas energy. 

Several numerical methods have been used to determine £. The Monte Carlo calculations of Carlson et al.[5] 
are generally regarded to be the most complete and gives £ = 0.44 ± 0.01. The same (or similar) results are 
showm in refs. [51 [7J [HI [H] some of which are £ = 4/9. Other authors report values of £ = 0.326 [TUJITT] while a 
recent result is £ = 0.25fT2]. 

It is a well-known fact that an interaction with large scattering length is separable. This paper is a report on 
results of calculations using a one term separable two-body interaction determined by inverse scattering from 
the phase-shift tv/2, i.e. with infinite scattering length and zero effective range. It is shown in Section 2 that 
the calculations are greatly simplified in the Unitary limit with an effective mass m* = 1. 

Numerical results are shown in Section 3. In the limit when the Pauli blocking Q — > 1, the theory reduces 
to the phase-shift approximation as shown in Section 4. A comparison with some of the results of other authors 
is shown in Section 5 and a short summary with comments is found in Section 6. 

2 Separable Interact ion; Formalism 

The use of a separable interaction in Nuclear Physics problems has a long history. It seems however that the 
first consistent calculation using inverse scattering techniques to construct a separable NN potential with an 
application to the nuclear matter problem was that reported in ref. [13j . A close agreement with calculations 
using the meson-theoretical potentials of Machleidt was found. Subsequent use was shown in ref. [TH [TS] 
relating to Vi ow -k etc. In the latter of these two last references the dispersion corrections and its relation to 
saturation was of primary interest. It is well known that for a two-particle system with a bound state at or 
close to zero energy the interaction can be represented by a separable potential. The method described below 
in which this separable potential is obtained by inverse scattering is therefore suitable when considering the 
problem at hand, large scattering lengths. One of the main problems in an inverse scattering calculation is the 
change in sign of the phasehift as a function of relative momentum. In the present case this is not an issue. 
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A rank 1 separable potential provides a sufficient and in fact precise description of the interaction in the 
Unitary limit. It is assumed to be attractive and given by 



V(k,p) = -v(k)v(p) (1) 
Inverse scattering then yields (e.g. ref [T3lfT6]) 

v 2 (k) = ^f s in5(k)\D(k 2 )\ (2) 
k 

where 



D(k 2 ) — exp 



2 k^ 

7T /„ fc 2 -fc' 2 



(3) 



where V denotes the principal value S(k) the phaseshift. A provides a cut-off in momentum-space. The effect 
of the cut-off will be exploited below. 
With S{k) = tt/2 one finds 

(A 2 - fc 2 )2 

and the interaction reduces to a constant for A 3> k, but — ► — oo for k — > A. 
The diagonal elements of the in-medium interaction is 

with 

wt-n-^jf^^P*-* (o) 

where P is the center of mass momentum, Q the angle-averaged Pauli-operator for pp-ladders (Brueckner 
approximation) and m* is the effective mass. One should note that the angle-averaging is exact in the effective 
mass approximation. 

The divergence of v 2 (k) for large k indicated after eq.(j4]) makes the numerical integration in eq.© somewhat 
complicated. This can be overcome by the substitution k' = Asin(t) to get 

1 r /2 m*Q(Asin(t),P)sin 2 tdt 
lG[k ' F) ~j2^ J (fc/A) 2 - smH (7) 

The effective interaction G(k,P) is in principle A-dependent. With m* = 1 and fc/A — > and Q — > 1 for 
A > kf, Ic{k,P) -> -1 + 0(1/ A). With w 2 (fc) — > 1/A one therefore finds G(k,P) independent of A for large 
A. If on the other hand one sets 1/A = in eq. |0 then G = K, where K is the reactance matrix as defined 
below. 

The (in)dependence of A is more clearly seen by using the following method applicable with m* = 1. One 
should note that both of the eqs ([6]) and ([7]) involve integrations up to A. A regularization that restricts the 
momentum-integration in eq. ([6]) to momenta < 2fc/ can be achieved as follows. 

With a separable interaction the Reactance-matrix K is defined by 

with 

w-^j^wh*** (9) 

where V refers to a principal value integration. Therefore K is real and its on-shell diagonal component is given 

by 

K(k) = -4TrtanS(k)/k. (10) 

In the Unitary case with S(k) — n/2, K(k) — > oo, implying Ix(k) — — 1. As stated above this is consistent 
with eq. (0 which for Q = 1 and m* = 1 also results in Iq = — 1 independent of k < A. 
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Eqs (0 and §E§ can be combined to get 



G(k,P) 



v 2 (k) 
I GK {k,P) 



(11) 



with 



I GK (k,P) = 



1 



(12) 



(2tt) 3 7 fc 2 - A:' 2 + i77 £an£(fe) 

The resulting form, ea. ljlip is very suitable for the problem at hand because in the unitary limit the last term 
in eq. (fl"2"l) — > 0. Furthermore, Q — ► 1 for fc' > 2fc/. The summation (integration) over intermediate states in 
eq. (fT2")l therefore only involves momenta < 2k f. One should observe that with Q defined for pp-ladders, Iqk 
is real for k < kf which is the case here. With hh-ladders one has pole-terms contributing to imaginary parts. 

Note that this regularization of the G-matrix equation is possible only because of the neglect of the dispersion 
correction i.e. with m* = 1. 

As shown by eq. v(k) is constant to any desired accuracy for momenta k < 2kf by choosing A large 
enough. Equation (TTTj) then simplifies. It is independent of the interaction because one can set v(k) = v(k'). 
The momentum-integration can then be done analytically as shown below. This simplification is not possible 
if using eq. ([6]) with integration up to the cut-off A where the interaction ([4]) diverges. It is also simpler than 
although in principle equivalent to eq. 0. 

As stated above, the Reactance matrix is real and so is the G-matrix (for occupied states) when defined 
with pp-ladders like in Brueckner theory. After dividing by v 2 (k') — v 2 k in eq. (I12p the integral is conveniently 
evaluated analytically One finds with a — and y — 



IcK(a,y) 



11 

IT 



1 + y + a * log 



1 + y-a 



1 + y + a 

To include hh-ladders there is an additional term 



+ - y 2 - a 2 )log 

2y 



(i + y) 2 



l-y 2 



+ kcotS(k). (13) 



hh(a,y) 



21 

7T 



2(1 - y 2 )* +a*log 



(1 — y 2 ) 2 + a 



(14) 



The diagonal G-matrix elements are 



G(a,y) = -4ir[I GK (a,y) - — + ^ofc 2 ]" 1 

(la A 



(15) 



where the two terms with scattering length a s and effective range vq drop out in the unitary limit but kept here 
for a discussion in Sect (5). 

With G(a, y) given, the potential energy per particle PE/A is 



3/c $ 



PE/A = —f 

k 2 Jo 



l-a 



8G(a iy )y 2 dy + - 



(l-a 2 )7 



4G(a,y)(l-y 2 - a?)ydy 



a da 



(16) 



The kinetic energy per particle, i.e. the uncorrelated fermi-gas energy is given by 

3 h 2 



E FG /A=^--k}. 

10 m 1 



The total energy is expressed in these units by 



Ej A = iE FG IA. 



2 This is a somewhat similar subtraction method as used to get eq. V(34) in ref.f |17|) and also used by S.A. Moszkowski in 
unpublished work [9]. 
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3 Numerical Results 



All results are for allowing the scattering length a s — > oo and effective range ro = 0. The integral I(a, y) = (i.e. 
G — » oo) along a line y = /(a) in the (a,y)-plane from approximately (0.84,0) to approximately (0.90,0.30). 
This complicates the numerical evaluation of the potential energy. This line is first extracted numerically by 
iteration at each meshpoint of the variable a to find /(a). The function I (a, y) is fitted to second order in y for 
each value of a in some interval Ay across this line f(a). The same techniques as used in standard principal 
value integrations is then used together with a shift in meshes so that one point will be located on the line /(a). 

The momentum-integration leading to the effective interaction was calculated analytically using eq. (|13p 
and also numerically from eq. ([7]) (with the same result). The integrations for PE/A, eq. (fT6|) . were done 
numerically. The result independent of density is 

£ = 0.24. 

It may be of some interest to see some results of approximations. So for example with 

9k t 

7(^)^1(0,0) = -^ 

one finds 

£ = 4/9 

which is a result also obtained by Steele [7J in the same approximation, but quite different from our £ = 0.24. 
If on the other hand one allows for a dependence on center of mass momentum by using 7(o, y) one finds 

£ = 0.515. 

One has to conclude that the momentum dependence on the effective interaction G(fc, P) cannot be ignored 
when calculating the energy of the many-body system. 

The above results are were obtained with a summation over pp-ladders as in Brueckner theory. One may 
inquiry as to the importance of hh-ladders for the present problem. By redefining the Q-operator by using 

Q = 1 — m — n-2 

instead of Q = (1 — nj.)(l — 712) {n being occupation-numbers) in the above equations one finds 

£ = 0.4 <-> 0.6. 

The uncertainty in the results is related to the pairing instability. Experience from nuclear matter calculations 
appears to be that this is to a large extent resolved by the Green's function method with integrations over the 
energy variable in the spectral functions. The above estimate of including the hh-ladders cannot be considered 
meaningful. A Green's function calculation is necessary in this case; results will be presented in a forthcoming 
report. 

To some approximation one might expect the pp- and hh-ladder contributions to be equal. This assumption 
was used in ref.[10j. Doing so here one finds 

£ = 0.56. 

But again, this result is not reliable. 

The convergence of G(k, P) when increasing A is not obvious from eqs ([5][7j. These eqs are however equivalent 
with eqs ()1 1I12[) which clearly do converge. That these two set of equations are indeed equivalent was also found 
numerically. But the last set of equations assume that to* = 1. It is found from numerical tests that eqs (07]) 
do not converge with increasing A so the problem with to* ^ 1 is not resolved here. As an example it is found 
that with m* = 0.9, £ -> 1 as A -> 10 3 fc/. 

4 Q — > 1; The Phase-shift approximation 

Two particles in a large box can be considered a limiting case of a many-body system. The effective interaction 
between two particles in a box having a relative momentum k is known to be given by [TH] [19l [20] [21} [22] 

G(k) = -An5(k)/k. 
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It is reasonable to expect a proper many-body theory of effective interactions to give this result in the limit 
Q — > 1. This limit has to be taken with caution. If one simply lets Q — > 1 in eq. (|12p and then takes the 
principal value one finds 

fcu 2 (fc) 
Gi ^" tan8(k) 

and 

G Q ^i(fc) = -Antan8{k)/k = K(k) 

where if (fc) was defined by cq. (fTO]) that relates to the /ree scattering of two particles rather than two particles 
in a box. (If instead inserting an irj in the denominator one obtains the complex T-matrix.) The correct limit 
is obtained after realising that when studying a many-body system one has to explicitly consider an enclosure 
of the particles in a large but finite box with a disrete rather than continuous spectrum. When extending the 
system so that sums can be replaced by integrations the transition from the discrete problem to the continuous 
has to be done with care. It has been shown [HI [19j [20] that when taking the limit Q — ► 1 one should in this 
case use 

- -> — + 7 



with 



Doing so in eq. (TT2"|) one finds correctly 



7 k( "d{k) tanS{k)^ 



G Q ^(k) =-4nS(k)/k. 



This may serve as a test of the equations. It was used in early work as an approximation to G{k) and referred to 
as the phase-shift approximation . It was believed to be a good approximation at low densities, substantiated by 
some numerical results. Of particular interest in relation to the Unitary problem is that it was used to calculate 
the binding energy for a neutron-gas at low density [231 124j . In this approximation of the effective two-body 
in-medium interaction one finds with 5{k) — -|, PE/A — — 4/3Efg giving 

quite different from any other result suggesting it to be a poor approximation here. The effect of the Pauli- 
blocking was however also calculated in ref.[23]. (The total energy of the neutron-gas reported there is some 
20% lower than reported in recent publications[5S].) Contrary to other beliefs the potential energy was found 
to be over-estimated by the phase-shift approximation by a factor varying between 1.75 — 1.95 for 0.1/m -1 > 
kf < 0.5/m -1 . These results were obtained with 1 5o phase-shifts and would not be directly applicable to the 
problem at hand. Using the same correction for the present Unitary problem would however give PE/A = 
(-0.76 <-> -0.68)^ G and 

£ = 0.24 <-> 0.32 

in fair agreement with our result £ = 0.24. 



5 Comparison with previous work 

The first publication relevant for comparing with the present work appears to be that of Baker [llj. He considered 
an attractive square-well potential with a radius c — > and an extrapolation of scattering length a — oc. The 
energy of the system was calculated in a pp-ladder approximation similar to the Brueckner G-matrix as used 
in the present report. It is however modified to avoid the Emery singularities [26]. The numerical evaluation of 
the energy using this 'R'-matrix gives a divergent result for /cpc — + present at all scattering lengths. A Pade 
approximant gave £ ~ 0.40. Baker also provides a series expansion of the ladder sum for c = 0. A [2/2] Pade 
approximant of this sum gives £ = 0.568 while the [1/1] gives £ = 0.326. 

Heiselberg [TU] has made extensive calculations both using a Galitskii resummation of hh- and pp-ladders 
resulting in £ = 0.33 and in a low order variational caculation resulting in £ = 0.46. The Galitskii method is 
somewhat similar to the present. One difference is that an average momentum was used when calculating the 
energy. As already mentioned above the hh- and pp-contribution were assumed to be equal which is another 
approximation. 
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The Monte Carlo calculations of Carlson et al[S] find a large pairing gap and a £ = 0.44 ± 0.01 includingthc 
pairing contributions. Without a pairing trial function (using a Slater determinant) they obtain £ = 0.54. □ 

Eq. (|13p for IcK{a,y) is the same function as f(n,s) (except for the ccrf-term) in the report by Steele, 
although with at least formally a very different method, using Effective Field Theory and Power Counting. 
[7] The expression for the potential energy (fT6| is consequently (in the limit of large sacttering length and 
with pp-ladders) also equal to that of Steele's. In his eq. (27) he lets f(n, s) — > 2 which is equivalent to 
our 7(0,0) = — - already considered as an approximation in Sect. 3. In this limit our results conequently 
agree giving £ = 4/9. This result is also obtained by Moszkowski[9 . Within the framework of Steele's work 
his approximation /(/t, s) — > 2 seems formally consistent with his expansion to order 1/2? . In the present 
work there is however no expansion other than in A and a s and the integrations over (a, y) gives a substantial 
correction from £ = 4/9 to £ = 0.24. 

Chen [8] like several other autors, finds £ = 4/9 but only in a low density expansion. He uses a relativstic 
approach motivated by the analogy with the infrared limit of Coulomb correlations. In the non-relativstic limit 
he finds the energy per particle to be a function of kf . 

Our result with pp-ladders, £ = 0.24 is appreciably smaller than in most reports. As mentioned in Sect. 4 it 
is to some extent substantiated by comparison with the neutron-gas calculations in ref. [24] . The calculations 
of Carlson et al 5J suggest that the pairing correction would give an even smaller £. The theoretical value closest 
to our result appears to be that of Lee [12] who calculated £ on the lattice with 22 particles in a periodic cube 
and found £ = 0.25(3). The estimate alluded to in Sect. 3 including the hh-ladders giving 

£ = 0.4 <-» 0.6 

lies in the range of most of the results but needs a better treatment by Green's function methods. 
There are experimental results reported between £ = 0.74 ± 0.07[2] and £ = 0.32to;io(2]- 



6 Summary and discussion 

This work addresses the problem of the energy of a zero temperature fermion gas in the unitary limit. The 
free-space two-body interaction is assumed to be separable. This assumption should in itself not affect the 
results that are expected to be independent of the shape or strength of the interaction defined only by the 
scattering length and effective range. The details of the interaction does not otherwise enter explicitly into 
the expression for the effective interaction or many-body energy. The separable interaction was in eqs. (1-4) 
determined by inverse scattering. Eq. ([4]) showed explicitly that in the unitary limit, a s — > oo and tq — > the 
interaction is constant, independent of momenta for k <C A. As A can be chosen arbitrarly large the interaction 
can then, with any chosen accuracy be considered a constant for all k < 2k f, 2k j being the maximum relative 
momentum in Q — 1 in eq. (|12j) . To perform the calculation of £ in the Unitary limit the effective interaction 
G can therefore be chosen to be independent of the interaction. 

The terms with a s and ro in eq. (|15|) do depend on the interaction. This equation would only be approxi- 
mate for finite values of these quantities and justified only to the degree that the interaction v(k) is constant. 
Numerical solution of eqs (1-3) shows that a necessary condition is that, as expected, tq = 0. 

The calculations are mainly analytical. The numerical integrations are simple. The result is as expected for 
a unitary limit independent of the assumed shape or strength of the interaction as well as of the density. The 
assumption of infinite scattering length sets the scale. 

The expression in eq. fl"6]) for the potential energy agrees with that of Steele using an EFT method and 
powercounting. This is not circumstantial, but ceratinly rooted in the fact that both methods, the separable 
potential and the EFT-power-counting rely on the nearly bound state with a pole near the real axis. There 
is a difference in final result in that the calculation here is carried a step further by doing the momentum 
integrations over occupied states. 

A selfconsistent inclusion of the hh- ladders and the related effect of spectral widths would be achieved in a 
Green's function calculation. In nuclear matter calculations these effects are found to be repulsive relative to 
the Brueckner results that include only pp-ladders with the potential energy increased by ~ 8%. [27] Applying 
such a correction here results in £ ~ 0.3. 

Green's function methods will be used in a forth-coming report on spectral functions and densities in 
momentum space. 

3 I am indebted to K.E. Schmidt to point this out to me. 
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There does not seem to be a consensus within the different theoretical calculations, ft would be expected 
that the most accurate are those of Carlson et al [5] . The result of the present investigation does not support 
those findings. 

The results within the present formalism are independent of density. This is expected to be general in 
the Unitary limit. At higher densities traditional many body theories predict higher order diagrams such as 
propagator modification by the mean field, included in standard Brueckner calculations to become important 
and. It is not clear how to deal with these problems in the Unitary limit. 

An example of such a problem is the effective mass. All calculations here are with to* = 1. The formalism 
used in the present calculations only works for this case. It is found however that the mean field is practically 
independent of momentum so that m* = 1 is realistic. 

BCS-pairing has not been included here. It is expected to lead to important corrections. Superfluid gaps 
have been calculated [TU] and are large for kpa > 1. QMC calculations show £ to decrease by about 0.1 when 
BCS-correlations are included in the trial wave- function [5]. 

The rather different results obtained for the coefficient £ in the Unitary limit both in experimental as well 
as theoretical reports suggests that this problem is still not resolved. 

It is somewhat intriguing that although £ is expected to be a universal constant its theoretical determination 
requires relatively complicated calculations. It is true that some estimates such as £ = 4/9 are the result of 
very simple assumptions and agree closely with the supposedly most accurate determination, ref [S]. Closer 
examination seem to suggest however (as in the present investigation) that these simple assumptions are not 
valid. 



It is a pleasure to thank Prof. Nai Kwong for many helpful discussions and Prof. Steve Moszkowski and 
Hcnning Heiselberg for some helpful suggestions and information. 
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